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A non-relativistic potential-model version of the factorization assump- 
tion, used in perturbative QCD calculations of hadronic form factors, is used, 
along with the Born approximation valid at high energies, to derive a remark- 
ably simple relationship between the impulse approximation contribution to 
the deuteron form factor at high momentum transfer and the high energy 
neutron-proton scattering amplitude. The relation states that the form fac- 
tor at a given value of is proportional to the scattering amplitude at a 
specific energy and scattering angle. This suggests that an accurate compu- 
tation of the form factors at large requires a simultaneous description of 
the phase-shifts at a related energy, a statement that seems reasonable re- 
gardless of any derivation. Our form factor-scattering amplitude relation is 
shown to be accurate for some examples. However, if the potential consists of 
a strong short distance repulsive term and a strong longer ranged attractive 
term, as typically occurs in many realistic potentials, the relation is found to 
be accurate only for ridiculously large values of Q. More general arguments, 
using only the Schroedinger equation, suggest a strong, but complicated, re- 
lationship between the form factor and scattering amplitude. Furthermore, 
the use of recently obtained soft potentials, along with an appropriate current 
operator, may allow calculations of form factors that are consistent with the 
necessary phase shifts. 
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I. INTRODUCTION 



The deuteron form factor A^Q"^) has been measured at Jefferson Laboratory and four 
momentum transfers up to Q'^ — 6 GeV^ [1], improving the SLAG measurements [2], and 
measurements at larger are planned. These efforts have caused much interest on improv- 
ing calculations of the form factors at higher values of up to about 11 GeV^. The best 
calculations are elegant in their use of the very latest realistic, high-precision nucleon-nucleon 
potentials. These potentials are based on using detailed knowledge of the long and medium 
range parts of the potentials and on using artful modeling of the short distance physics. 
Typically, the parameters of the potentials are tuned to obtain an accurate reproduction 
of the measured phase shifts up to 300 McV laboratory kinetic energy, T. Increasing the 
range of energies of the validity of the potential should increase the ability of the potential 
to describe those aspects of the deuteron wave function which enter at high momentum 
transfer [3]. 

Indeed it seems reasonable to expect that describing a form factor at a given momentum 
transfer, Q^, would require a reproduction of the large angle NN scattering amplitude at 

where uin is the nucleon mass. Qualitatively this is because the relative momentum, which 
dominates the overlap integral for the form factor is ~ Q/2. The implication of (1) is that 
for > 1 GeV^ one needs to use a potential that successfully describes the neutron-proton 
phase shifts for T > 500 MeV. To our knowledge, no such consistency check of the potential 
used to compute the deuteron wave function has been made. 

One might think that a simple kinematic relationship such as Eq. (1) might not apply 
because the nucleons in the deuteron are bound. Therefore we need to demonstrate explicitly 
that there is a strong relationship between the form factor and the scattering amplitude at 
large energies and large angles. Indeed, the specific result (see Eqs. (15) (17) below) 
differs quantitatively from Eq. (1). This is obtained using the simplest possible dynamics: 
non-relativistic spin-less nucleons interacting with an energy-independent local potential. 
Glearly, it is not our purpose to be realistic. Instead we merely wish to point out that, 
under certain conditions, the form factor can be proportional to the scattering amplitude. 

The derivation of the relation between the form factor and the scattering amplitude pro- 
ceeds in Sec. II by applying the factorization approximation commonly used in perturbative 
QGD derivations of hadronic form factors along with the first Born approximation expected 
to be valid at high energies. The requirements for each approximation are investigated and 
correction terms obtained. The accuracy of the approximations and the resulting form fac- 
tor scattering amplitude relation are studied using simple interactions: attractive Goulomb 
potential and sum of attractive and repulsive square wells in Sect. Ill, and some implications 
for other models are discussed. A detailed numerical study using the Malfliet-Tjon potential 
is made in Sect. IV. The specific approximations used in Sect. II appear to be marginally 
successful, if the potential obtains a total weak attraction by combining strong repulsion 
at short distances with strong attraction at larger distances. Therefore, a more general ar- 
gument, using only the Schroedinger equation, is presented in Sect. V. A discussion of the 
implications of our results is presented in Sect. VI. 
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The reader may immediately question the use of non-relativistic dynamics, because many 
specific relativistic effects are known. For a recent review, see Ref. [4]. However, such 
dynamics are not irrelevant at T = 500 MeV. Furthermore, the use of relativistic light 
front dynamics shows that relativistic dynamics is not very different from non-relativistic 
dynamics: two-nucleons dominate, there is a wave equation, and the specific relativistic 
effects in the deuteron are not very large unless is very high. The specific differences 
between the non-relativistic and light-front approaches are relatively well-understood and 
lead to a small easing of the constraint (1) for > m^. Therefore, we turn to the necessary 
derivation without further apology. 



II. BASIC IDEA 

The deuteron wave function ip is defined by the Schroedinger equation: 

where /j, = /2. The form factor is given in terms of the momentum-space wave function 
as 

F{q) = j d^pij{p)ij{\p + q/2\)F{q) = J d^pij{p)ij{\p + q/2|). (3) 

In the widely-employed Breit frame = Q^. For large enough values of q, and for a 
potential (expressed in momentum space) that which decreases as a power of q, the integral 
may be simplified because q can be much larger than the typical values of p for which the 
wave function is near its largest value. The deuteron wave function is known to have a 
limited momentum content, being essentially for momenta greater than about 600 MeV/c. 

We aim for a simplification of Eq. (3) that is vahd at large momentum transfer. The 
basic assumption is that, for large enough values of q, one can regard q ^ p even though 
the integral extends over all values of p. This leads to the factorization approximation 
commonly used in perturbative QCD calculations of hadronic form factors. A pedagogic 
discussion of the technique we employ is presented in Ref. [5]. There are two regions for 
which the integrand of Eq. (3) is largest: p for which the first wave function is large, 
and p — q/2 for which the second wave function is large. The regions contribute equally 
so that we may say 

FiQ') ^ FaiQ') = 2^(g/2) J £p (4) 

with the integral over all momentum being proportional to the coordinate space wave func- 
tion evaluated at the origin: 

^{r = 0) = / ^{p). (5) 



One may also derive (4) using coordinate space arguments and determine the leading 
correction term. By regarding p <^ q/2 one finds [6] the approximate form factor and its 
leading correction term: 
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d^pi^ip) J j—^ij{r)exp{i{p + q/2)-r)) 
d^r 



I d^p^ip) I -^V'(r)(l-pV/6)exp(iq/2-r) 



^ Fa{Q') - {2nf/'V'^^;{r = O)S^(^). (6) 

If the integrals exist, and if the wave functions fall as a power of q, the ratio of the second 
to the first terms of Eq. (6) is proportional to ^ j^I^p^^^-^ I which vanishes for sufficiently 
large values of g^. 

According to Eq. (4), the momentum transfer dependence of the form factor is obtained 
from the wave function ip{q/2), for sufficiently large values of Q"^. Thus we attempt to obtain 
^^{q) from the momentum space version of Eq. (2): 

i;{q/2) = f d'p (q/2|y|p) ^(p). (7) 

We again use the idea that q can be much greater than p, so that the wave function 
can be approximated as 



i:{q/2) ^ ^(q/2|V|0) / d^pilj{p) = -^{q/2\V\0){27r)h{r = 0). (8) 
g"^ J q'^ 



This approximation seems very natural, if one is used to the ideas of perturbative QCD. In 
particular, the entire momentum transfer is taken up by a single action of the potential, so 
that the important positions in coordinate space region are those for which the potential 

has its greatest variation. 

It is useful to examine Eq. (8) from the view of coordinate space. One has 

j rf-V(q/2|V|p) ^(p) = j j^^exp{-iq/2 ■ r))y(r)V^(r). (9) 

One may now observe that the approximation Eq. (8) relies on replacing the product 
V{r)ilj{r) by V{r)ip{0). This is allowed only if the potential has a much more significant 
variation than the wave function for positions near the origin. 

The expression Eq. (8) can be used to obtain the form factor from Eq. (3). The result is 

F,{Q') ^ WV\0) $^ (10) 



$ = 



= jd^pi^{p). (11) 



The approximation (10) has been verified numerically in Ref. [7] for the case of the semi- 
realistic Hulthen model. 

The content of Eq. (10) is that the form factor is proportional to the square of the bound 
state wave function at the origin times the Fourier transform of the potential: 
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{ci/2\V\0) = d'r e'^^-V(r) (12) 



f r, sin igr , , , ^ 

j r^dr 1 y y(r). (13) 



(27r)3 J 

The desired relation between the form factor and the scattering amphtude can be ob- 
tained by reahzing that at energies E much greater than the characteristic strength of the 
potential, the first Born approximation is valid and the scattering amplitude for a cm angle 
6*, /e(6') = /^(k, k'), E = A;^/2/i = /c'^/2/x, k • k' = k'^cosO is also a Fourier transform of 
the potential. The validity of the first Born approximation, at high energies, for energy- 
dependent local potentials, is well-verified in many quantum mechanics textbooks. In that 
case, 

fE{e) ^ / d're^^^-'^'y'^Vir). (14) 

The relation between /e(6') and (q/2|V|0) is obtained by specifying 

^q'^ = 2A;2(1 - cos^) = AiiE{l - cos 9) = Tmjv(l - cos^). (15) 

Imposing the relation (15) immediately yields 

fE{9) = -47rV(q/2|l^|0), (16) 

and with Eq. (10) 

P(Q') = :^2Ue) ^\ (17) 

This is the result we have been seeking. The key point is that the form factor is expressed 
as the high-energy scattering amplitude times well-defined factors, implying that a correct 
calculation for the form factor can only be achieved in models in which the scattering 
amplitude is accurately reproduced. 



III. TOY MODELS 

The arguments used in the previous section are analogous to those presented in deriva- 
tions of QCD factorization theorems. However, these have not been used often for nuclear 
targets. We therefore discuss two simple examples. 



A. Coulomb binding 

We take the potential to be 



V{t) = (18) 
r 
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with the exact wave function given by 



with a — l/{^g^). The form factor F{q) of Eq. (3) is given by 



(19) 
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(20) 



{a?q^ + 16)^ 



The wave function in momentum space is 



V'(p) 



47riV 
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(21) 



(27r)3/2 




2 ' 



a 



so that the approximation of Eq. (4) yields 



Fa is close to F for q^a^ ^ 3. Using = a^/2^ and the deutcron binding energy and 
H = Mjq/2 means that all that is required is g ^ 75 MeV/c. The approximation of Eq. (10) 
yields 



which valid under similar conditions. 

The next task is to determine the conditions that the first Born approximation Eq. (14) 
be valid. It is well known that, for the Coulomb interaction, this approximation reproduces 
the exact scattering cross section. However, the correct scattering amplitude is complex, 
while the approximation Eq. (14) gives a real result. Thus the condition that Eq. (14) be 
valid is the condition that the s-wave scattering phase shift be small. This is the condition 
that g^MN/{2k) < 1, or A;a > 1. Since a = l/(45MeV/c, one needs only k » AbMeV/c. 

Thus the approximations of the previous section are easily satisfied for dynamics defined 
by the Coulomb potential. More generally, it is reasonable to expect that, if the potential is 
local, purely attractive with a significant gradient at the origin and magnitude determined 
by the very small dcuteron binding energy (2.2 McV), the approximations needed to reach 
the form factor scattering amplitude relation (17) are well satisfied. In such cases, the 
potential varies more rapidly than the wave function for positions near the origin, and the 
potential is weak enough so that the first Born approximation can be accurate for reasonable 
momentum transfers. 

B. Short- ranged repulsive square well plus long ranged attractive square well 

The nuclear force is repulsive at short distances and attractive at long distances between 
nucleons. Thus we consider the model defined by 
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(23) 




V{r) = Voe{E^ - r) - Vie{R, - r)e{r - Ro), V^o,i > 0. 



(24) 
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We take i?o = 0.4 fm and Ri = 1.5 fm. Then the observed deuteron binding energy is 
reproduced using Vq = 0.302 fm^^, Vi = 29.18 fm~^ = 5.757 GeV. This corresponds to 
a hard core repulsion. The wave function t/j^r) is given (up to an overall normalization 
constant) by 



, r<Ro 



sihRRq sinh7r 
sin KRi sinh 'jRo ' 

sin Kr „ „ 
-, Ro<r < Ri 



sin KRi 

exp( — /•/«) 
exp{—Ri/a) 



,r>Ri 



(25) 



(0 
1^ 



-1 




6 7 8 

Log(Q/fm-^) 



FIG. 1. Ratio of approximate to exact form factor. 

The exact {F{q)) and approximate {Fa{q) Eq. (4)) form factors are compared by dis- 
playing their ratio in Fig. 1. The are some wild fluctuations arising from nodes of F{Q), 
but these are not the most interesting feature. The true difficulty is that the approximation 
becomes valid only for extremely large values of q > 100 fm~^ . Furthermore, the ratio 
Fa{q)/F{q) is =0.87 instead of unity. That difficulties exist can be seen immediately by 
comparing Eqs. (4) and (25). The coordinate-space wave function vanishes at the origin, 
so the leading approximation vanishes! This signals a breakdown of the derivation. The 
point is that the factorization Eq. (4) depends on the potential varying more rapidly near 
the origin than the wave function 

One might then be amazed that the approximate form factor is even close to the exact 
one. That this occurs can be understood from the discussion below Eq. (6). The action of 

does not introduce an extra factor of because the momentum space wave function 
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contains terms ~ — ^2 . 

This section contains two simple models. The approximations work magnificently if the 
Coulomb potential is used, but fail miserably with the nuclcon-forcc-motivated model of 
Eq. (24). Let's examine a more realistic model to see if our approximations are relevant for 
understanding deuteron physics. 



IV. MALFLIET-TJON POTENTIAL 

The short-distance repulsion and longer-ranged attraction are simulated via potentials 
of the Yukawa form in the Malfliet-Tjon potential [9]: 

V{r) = -Xa + Xr . (26) 

r r 

The parameters were chosen to reproduce the deuteron binding energy, scattering length, 
effective range and phase shifts up to a laboratory energy of 300 MeV. This potential is 
much smoother than the square wells of Eq. (24). 

Our aim is to study the approximations given in Eqs. (4,8) and the Born approximation 
of Eq. (14). The comparison between the exact form factor F{Q'^) and the approximation 
F„(g2) of Eq. (4) is displayed in Fig. 2. 




log(Q fm) 



FIG. 2. Ratio of approximate to exact form factor for the Malfliet-Tjon deuteron wave function. 
The solid curve compares the approximation of Eq. (17) with the exact form factor. The dashed 
curve compares the approximation of Eq. (4) with the exact form factor. 

The approximation does work, but only for huge values of Q. The value of Q ~ 140 
fm~^, or a; 900 (GeV/c)^ for its natural logarithm to reach the value of 5. One may 
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search for the cause of this bizarre hmit by examining the correction term shown of Eq. (6). 
This is governed by the kinetic energy operator acting on the wave function at the origin. 
The Schroedinger equation that this essentiaUy the potential times the wave function, with 
a product that varies as {Xr — A^) near the origin. The strong nature of the repulsive term 
Xr — 7.41 causes the expansion shown in Eq. (6)to converge very slowly. 

One may gain further insight by studying the relationship between the exact and ap- 
proximate wave function of Eq. (8). As shown in Fig. 3, the approximation attains va- 
lidity only at supremely large values of the relative momentum. Here accuracy requires 
V{r)ijj{r) ^ y(r)^(0) for values of r near the origin. For the Malfiiet-Tjon potential 
V{r) ~ 1/r and ip{r) approaches 0, in contrast to the Coulomb wave function. From this, 
the failure of the approximation seems natural. Moreover, many realistic potentials behave 
in a similar fashion. 



1 — \ — \ — r 




2 4 6 8 10 

log(p fm) 

FIG. 3. Ratio of approximate to exact wave function for the Malfliet-Tjon deuteron wave func- 
tion. 

The final requirement needed to achieve the result (17) is the vahdity of the Born approx- 
imation. Our comparison of the exact and Born approximation to the forward scattering 
amplitude is shown in Fig. 4. The partial wave expansion is used unless the energy is large 
enough (lab energy greater than one or two GeV) for the eikonal approximation to repro- 
duce the exact scattering amplitude. Fig. 4 shows that the Born approximation does become 
valid, but only at absurdly high energies. This figure shows results only for forward scat- 
tering, whereas our relation (17) requires high momentum transfer. The convergence to the 
Born approximation occurs for higher energies for large scattering angles, if the Malfiiet-Tjon 
potential is used. 

The slow approach to the Born approximation can be understood from the eikonal for- 
malism in which the scattering amplitude is expressed in terms of an integral 
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fOO 

fM^-ik hdh M2khs\u(d/2))(e'^^^^ -I), (27) 
Jo 

with the phase shift function x(6) given as 

x{h) = dzV{V¥T^). (28) 

For the Malfliet-Tjon potential, x(6) is given by 

Til 

X{b) = jiXuMbi^n) - XAMbi^A)), (29) 

which is small only for very large values of the relative momentum k. While the overall 
potential is able to produce only one weakly bound state, the weak attraction arises from 
the cancellation of two very strong terms of opposite sign that generally causes large values 
oixib). 




FIG. 4. Ratio of approximate to the real part of the exact forward scattering amplitude for the 
Malfliet-Tjon potential. E is a laboratory energy. 

The factorization approximations for the Eqs. (4,8) and the Born approximation of 
Eq. (14) work for sufficiently large values of the transferred momenta. A reasonable reader 
might examine the figures 2-4 and conclude that "sufficiently large" means too large to be 
relevant for real experiments. However, success at asymptotically high momenta indicates 
that the key assumption Eqs. (4,8) that the form factor is proportional to the wave function 
at the origin is a very stringent requirement that is especially difficult to satisfy for a wave 
function that vanishes at the origin. Furthermore, high momentum transfer scattering may 
proceed via a set of small momentum transfer processes under certain circumstances. Thus 
the condition that the Born approximation (scattering by a single action of the potential) 
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be valid is also a very strong one, as is spectacularly manifest in Fig. 4. Truly magnificently 
high energies are required for the Born approximation to be valid. 

While the relation Eq. (17) suggests a close connection between the form factor and the 
scattering amplitude it is worthwhile to ask if one can derive a connection between the form 
factor and the scattering amplitude that does not directly involve the wave function at the 
origin and single-scattering assumptions. The desired connection is displayed in the next 
section. 



V. GENERAL APPROACH 

The central idea is that the wave function is dominated by its low momentum components 
and that high momentum components are obtained by at most one high momentum transfer 
operation. Wc intend to manipulate the Schroedinger equation to derive an explicit relation 
between the form factor and the scattering amplitude. Projection operators are used to 
develop the necessary formalism. Let 

g = /dV(A-p)|p)(p|, p = Jd'pe{p-A)\p){pi (30) 

with P a projection operator on high momentum states and A a parameter, approximately 
600 MeV/c or higher, denoting the separation between the high and low momentum transfer 
regions. Then the Schroedinger equation can be expressed in terms of low, Q\ip) = \'4')q, 
and high, P\ip) = \4')p, momentum components: 

^ {Vpp\^)p + Vpq\^)q) (31) 



—sb— Tpp 

Wq = {VqqWq + VqpWp) , (32) 

-€b - -LQQ 

where the notation PTP = Tpp, PV P = Vpp, etc. is used for the kinetic energy T and 
potential energy operators. Substitute Eq. (32) into Eq. (31) to obtain: 

= (vpp + VpQ ^ T^^Q^] 1^)^- (33) 

-cb - -ipp \ -eB - iQQ- Vqq J 

Consider the complete eigenstates of the Hermitian, energy-independent operator Tqq -\- 

Ek\<j)k)^-^ = {TQQ + VQQ)\<Pk)^-^ (34) 
E^\c^n) = iTQQ + VQQ)\c^n) (35) 

in which \(j)k)^~^ is a scattering state with incoming boundary conditions, Ek = k'^/2fi, and 
En 7^ —eB, \4>n) are the energies and wave functions of any bound states that may exist. The 
use of completeness in Eq. (33) leads to 

— efi — ipp V J —eB—Ek n —^B—En I 
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The form factor depends on the high momentum components of the deuteron wave 
function, so consider the quantity (q/2|-?/')p for Qjl > A. Then {ci/2\ijj)p = (q/2|?/'). 
The second term on the right-hand-side of Eq. (36) is determined by the matrix elements 
{q/2\V\(pk)''~\ which are scattering amphtudes for off-energy-shell kinematics [18]. These 
matrix elements are the terms which contain the relation between the form factor and the 
scattering amplitude. To see this, take the parameter A to be just a bit less than q/2. 
The largest of the matrix elements (q/2|y|0fe)*^^^ will be those for which k is just a bit 
less than A. In this case, the kinematics are nearly on-shell. Transition matrix elements 
are continuous, so that it is reasonable to expect that reproducing the observed scattering 
amplitudes at high energy is important for reproducing the important features of the high 
momentum components of bound state wave functions. 

VI. DISCUSSION 

The use of two factorization approximations, Eqs. (4) and (8) combined with the first 
Born approximation, Eq. (14), lead to a statement (17) that the form factor is proportional 
to the scattering amplitude times the square of the coordinate space wave function at the 
origin divided by Q^. A relation very similar to Eq. (17) was derived long ago [8] using 
a scale-invariant, six-quark model. It is also true that the relation between the deuteron 
form factor and the scattering amplitude has also been the subject of Ref. [10] in which 
dispersion relations are used to compute the deuteron charge form factor with experimental 
phase shifts as the essential input. Thus one sees the close relation between F(Q^) and 
Je from a variety of different approaches: non-relativistic dynamics, light-front dynamics, 
quark models and dispersion relations. 

But the validity of Eqn. (17) depends on the use of factorization approximations and the 
Born approximation that might not be valid for realistic nucleon-nucleon potentials. The 
analysis of toy models. Sects. II, and the Malfliet-Tjon potential Sect. Ill, suggests that only 
a soft potential can be expected to satisfy the conditions necessary to obtain the result (17). 
Until recently, the idea that a soft potential could also be realistic was only a faint hope. 
But recent work, using effective field theory have introduced a set of soft-realistic potentials 
[11], [12], [13]. It is possible that these potentials, which do not contain those features of the 
Malfliet-Tjon potential that violate the necessary approximations, could satisfy Eqn. (17). 

It is worthwhile to suppose that Eq. (17) could be accurate for some soft realistic poten- 
tial. Then we may estimate the kinetic energies required for reproducing the phase shifts. 
Since we are dealing with the S-wave deuteron wave function we can take cos^^ = 0, to 
correspond to the maximal momentum transfer for identical particles. In this case, Eq.(15) 
gives 

= AtunT. (37) 

The use of this is simple. Suppose the phase shifts are well described up to about T — 350 
MeV, which is the standard stated upper limit. Then one can calculate the form factor up to 
only Q'^ ^ lA GeV^. But modern high-precision NN potentials actually describe the data 
up to laboratory energies of about 1 GeV [14,15]. Thus one may compute the form factors 
up to ~ 4 GeV^. However, present measurements reach = 6 GeV^, and there are 
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plans to reach higher values. Thus the constraints we present may be relevant for present 
and future measurements. 

An additional worry, is that the analysis of the previous paragraph might require a 
relativistic treatment. If we instead were to use hght cone models as in [16,17] we would 
find that up to ~ the light cone fractions of the nucleons are approximately equal 
and the relation Eq. (15) holds, but with T = 2E + E'^/2mN- At larger values of Q^, the 
increase of the effective invariant energy with decreases somewhat. In any case, one 
can see that for > 2 GeV^ one reaches the region where masses in the intermediate 
state exceed 3 GeV and the legitimacy of the two nucleon approximation becomes highly 
questionable, as discussed in Ref. [16]. Then one would need to include either new hadronic 
degrees of freedom in the dcuteron wave function: ttNN, AA, or to account explicitly for 
quark-gluon degrees of freedom. 

The implications of our result (17) is that potentials used to obtain the deuteron wave 
function should be tested by computing the corresponding phase shifts. If one wants an 
accurate calculation, the phase shifts need to be correctly obtained up to kinetic energies 
given by Eq. (37). This requirement appears, according to the arguments of Sect. V, to be 
more general (see Eq. (36) than accuracy of the asymptotic relation between the scattering 
amplitude and the form factor, (17) as the slow onset of asymptotia occurs, in part, because 
the wave function at the origin is severely suppressed by strong short distance repulsion, 
and because single hard scattering is not strongly dominant over the possibility of achieving 
a high momentum transfer as a result of several small momentum transfer processes, large 
extent due to cancellations between the contribution dominating in the asymptotic region 
and other contribution. All together this suggests that potentials currently employed to 
compute deuteron form factors be tested for consistency with high energy phase shifts. A 
soft potential [11], [12], [13] and an appropriately derived current operator should also be 
used. 
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